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This article introduces a dynamical Ginzburg-Landau phase transition/separation model for the 
mixture of liquid helium-3 and helium-4, using a unified dynamical Ginzburg-Landau model for 
equilibrium phase transitions. The analysis of this model leads to three critical length scales L\ < 
L-2 < L3, detailed theoretical phase diagrams and transition properties with different length scales 
of the container. 

PACS numbers: 



I. INTRODUCTION 

Landau's legacy on phase transition has been a singular driving force for our recent work on developing a general 
dynamic transition theory and a unified approach for both equilibrium and non- equilibrium dynamic phase transitions. 
We hope that the work in this paper is in the spirit of Landau's legacy. 

Superfluidity is a phase of matter in which "unusual" effects are observed when liquids, typically of helium-4 or 
helium-3, overcome friction by surface interaction when at a stage, known as the "lambda point" for helium-4, at 
which the liquid's viscosity becomes zero. 

The main objectives of this article is to study A-phase transitions of liquid helium-4 and and phase separations 
between liquid helium-3 and liquid helium-4 from both the modeling and analysis points of view. 

In the late 1930s, Landau proposed a mean field theory of continuous phase transitions. With the successful applica- 
tion of the Ginzburg-Landau theory to superconductivity, it is nature to transfer something similar to the superfluidity 
case, as the superfluid transitions in liquid 3 Hc and 4 Hc arc of similar quantum origin as superconductivity. Unfortu- 
nately, we know that the classical Ginzburg-Landau free energy is poorly applicable to liquid helium in a quantitative 
sense, as described in by Ginzburg in pQ. 

As an attempt for this challenge problem, we introduces a dynamical Ginzburg-Landau phase transition/separation 
model for the mixture of liquid helium-3 and helium-4. In this model, we use an order parameter ip for the phase 
transition of liquid 4 He between the normal and superfluid states, and the mol fraction u for liquid 3 He. As u is a 
conserved quantity, a Cahn-Hilliard type equation is needed for it, and a Ginzburg-Landau type equation is needed 
for the order parameter. The interactions of this quantities are built into the system naturally by using a unified 
dynamical Ginzburg-Landau model for equilibrium phase transitions, where the dynamic model is derived as a a 
gradient-type flow as outlined in the appendix. 

This analysis of the model established enables us to give a detailed study on the A-phase transition and the phase 
separation between liquid 3 He and 4 He. In particular, we derived three critical length scales L\ < L2 < L3 and 
the corresponding A-transition and phase separation diagrams. The derive theoretical phase diagrams based on our 
analysis agree with classical phase diagram, as shown e.g. in Reichl [5] and Onuki [TJ, and it is hoped that the 
study here will lead to a better understanding of mature of superfluids. Finally, we remark that the order of second 
transition is mathematically more challenging, and will be reported elsewhere. 

One important new ingredient for the analysis is a new dynamic transition theory developed recently by the 
authors [2J [3J. With this theory, we derive a new dynamic phase transition classification scheme, which classifies 
phase transitions into three categories: Type-I, Type-II and Type-Ill, corresponding respectively to the continuous, 
the jump and mixed transitions in the dynamic transition theory. 
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II. MODEL FOR LIQUID MIXTURE OF 3 HE - 4 HE 



Liquid 3 He and 4 He can be dissolved into each other. When 3 He -atoms are dissolved in liquid 4 He and the density 
of 3 He increases, the A-transition temperature T\ decreases; see the liquid mixture phase diagram of 3 He - 4 He (Figure 
[lj , where X — 77,3/(77,3 + 714) , 713 and 72.4 are the atom numbers of 3 He and 4 He respectively. When X = 0, T\ — 2.17K, 
where A-phase transition takes place and liquid 4 Hc undergoes a transition to superfluid phase from the normal liquid 
phase. When X = 0.67 and temperature decreases to T = 0.87K, i.e., at the triple point C in Figure [l] the liquid 
mixture of 3 He - 4 He has a phase separation. 
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FIG. 1: Liquid mixture phase diagram of 3 He - 4 He . 

Let the complex valued function t/j = tpi + iip2 describe the superfluidity of 4 He, and u be the density of 3 He, which 
is conserved, i.e., 



/ udx = c (c > is a fixed number). 



The Ginzburg-Landau (Gibbs) free energy is taken in the following form 

-k x h 2 



G(ip,p ni u) 



2m 1 1 1 2 ITI 4 111 2' 1 2 3 4 

where h is the Planck constant, m is the mass of helium-4 atom, and the coefficients satisfy 

^1,^2,72,^2,^3,^4 > 0. 



dx. 



(1) 



(2) 



(3) 



When -0 = 0, the liquid mixture of helium-3 and hclium-4 is a binary system, and G(0, u) stands for the Cahn-Hilliard 
free energy. By the standard model ( A9 1, from ([!]) and ^ , the equations governing liquid mixture of 3 He - 4 He are 
given as follows: 



dt m 
(hi 

di 



(4) 



-k 2 A 2 u + A[v x u - i/ 2 u 2 + v 3 u 3 + Vili^l 2 ] ■ 
These equations have a physically sound constant steady state solution given by 

(V,U) = (0,TZ°), 

where u° > is the density of helium-3 in a homogeneous state. For simplicity, we assume the total density p = 1. 
Then the system control parameter X, the mol fraction of 3 He, becomes 

X = u a with < X < 1. 

Now we consider the derivations from this basic state: 
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and we derive the following equations (drop the primes): 

fiiAip — Xiip — 72|?A| 2 V' — 2v4Uijj, 



at 

— = -n 2 A 2 u + \ 2 Au + A[(-v 2 + 3v 3 X)u 2 + v 3 u 3 + u 4 \il>\ 2 ], (5) 
udx = 0, 



with the Neumann boundary condition 



where mui = k\h 2 /m, [i2 = k%, and 



d 

(u, Au,tp, p n ) = on c?Sl, (6) 



dn 



(7) 



Ai = 71 + 2v 4 X, 
\ 2 = v 1 - 2v 2 X + 3v 3 X 2 . 

III. PHASE TRANSITIONS OF 3 HE - 4 HE MIXTURES 

We now apply equations Q to study the phase transitions of 3 He - 4 He liquid mixtures. 

A. Critical parameter curves 

For simplicity we only consider the special case where the container is a rectangle: 

= (0, L) x (0, 1) 2 C M 3 for L > I, 

and the control parameters are the temperature T and the mol fraction X, and the length scale L of the container f2. 

Physically, by the Hildebrand theory (see [H]), in the lower temperature region and at p = 1 atm, the critical 
parameter curve X 2 = is equivalent to 

T= 2 *(l-X)X-a , (8) 

where R is the molar gas constant and a > is a constant. Here <7o > is small correction term. The original 
Hildebrand theory leads to the case where (To = 0. However, as we can see from the classical phase separation of a 
binary system, the Hildebrand theory fails when the molar fraction is near or 1, and the correction term added here 
agrees with the experimental phase diagram as shown e.g. in Figure 4.13 in Reichl [8 . 
Furthermore, by Q and ([8]), we have 

p 1 ~0 1 {T + a o ), 3^3-2^ 2 , (9) 

Consider eigenvalue problem of the linear operator in (|5| : 

ViA - Ai 

Hi A — Ai I I V>2 I = P I ife I • (10) 



-A(/i 2 A - A 




Here ip = ipi + ip 2 . It is known that the first eigenvalue and eigenvector of the Laplacian operator with the Neumann 
condition and zero average are A = ir 2 /L 2 and u = costtxi/ L. Thus, the first two eigenvalues and their corresponding 



eigenvectors of (10) are given by 

f/3i = -Ai = -( 7 i + 2z/ 4 X), 



(ipi,ip 2 ,u) = (1,0,0) and (0,1,0), (U) 

*~p0£ + *)~£0£ + «- 2 -*+ , 12 , 

(lpl,1p2,U) — (0, 0, COSTTXl/L), 
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As for v\ , the parameter 71 is approximately a linear function of T. Phcnomenologically, we take 

7i = -<t 1 +# 2 T (ff 2 ,ffi>0). 
Then, by ^ and (ll|-(13l, the critical parameter curves in the TJT-plane are as follows: 

(Ti — 2v±X <j\ 



(13) 



h ■ Pi = 
Z 2 : A = 



Tel 



6, 



(i-x), 



c2 R 1 ; 6>iL 2 



(14) 



Here we have assumed that a\ — 2v±. This assumption is true at X = 1 as required by T c i = at X = 1, and is true 
as well at other values of X as we approximately take o~\ and as constants. 
The critical parameter curve l\ (Pi — 0) is as shown in Figure |2ja). Let 
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then when L > L the critical parameter curve I2 (P2 — 0) is as shown in Figure |2 



(15) 





(a) (b) 

FIG. 2: (a) /3i < in region Ai, and f3i > in A 2 , (b) /? 2 < in region Si, and /3 2 > in B 2 . 



B. Transition theorems 



The following is the transition theorem for liquid mixture 3 He - 4 He. Fro this purpose, we introduce a few length 
scales as follows: 
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(17) 
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> cr 



(18) 



otherwise. 
We remark here again that o~q is small as a correction term in the Hildebrand theory as mentioned earlier in (T8J) 



Theorem 1 Let T c \ and T C 2 be given in (14). For equations |5}) ; we have the following assertions: 
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(1) When L < L\ the system has only the superfluid phase transition (i.e., the X-phase transition) at T = and 
the TX -phase diagram is as shown in Figure^a). 

(2) Let <7i/#2 < 4a/ R. Then L\ < L2, and if L\ < L < L2, then there are two numbers < X\ < | < X2 < 1 
given by 



X\.2 — 



1 



V a V^i72 O1L 2 



(19) 



such that the following hold true: 

(a) if < X < Xi or X 2 < X < 1, the system has only the X-phase transition at T = T\; 

(b) if Xi < X < X2, then the system has the X-phase transition at T = T c \, and has the phase separation at 



T* = T C 2 + < T c \. Moreover, the TX -phase diagram is as shown in Figure |£| 

(3) Let <7i/#2 < 4a/i?. For any L2 < L < L%, the phase diagram is as shown in Figure [7J For any L3 < L, the 
phase diagram is as shown in Figure [5[ 

(4) Let ai/02 > ia/R. Then for any L > L\, Assertion (2) holds true. 

(5) The X-phase transition at T = T]. is Type-L, which corresponds to second-order transition. 




FIG. 3: Region D\ is normal 3 He - 4 He , D2 is the superfluid phase, and D3 is the phase separation region. 




FIG. 4: Region D\ is normal 3 He - 4 He , D2 is the superfluid phase, and D3 is the phase separation region. 

It is worth mentioning that the order of second transitions crossing f3i = is an interesting problem, which will be 
analyzed in a forthcoming paper. 

Proof of Theorem [TJ We proceed in several steps as follows. 

Step 1 . It is easy to see that the space 



E={(u,^) = (0,tJj)\ VeC} 



OX, X" x 2 1 X 

FIG. 5: Region D\ is normal 3 He - 4 He , D2 is the superfluid phase, and D3 is the phase separation region. 



is invariant for ([5]). Therefore, the transition solutions of Q from the critical parameter curve l\ (Ai = 0) must be in 
E, which corresponds to the superfluid transition. On the other hand, restricted on E, equations Q are equivalent 
to the following ordinary differential equation: 



dtp 
~dt 



-Ai^-tsIVIV (20) 
It is then easy to see that Assertion (5) holds true, and the bifurcated solutions consist of a circle given by 



I [0,2tt]}. 

Step 2. We now consider the second transition. For this purpose, let the transition solution of Q from T = 
(i.e., from Ai = 0) be given by 



(0,^(T)) = (0, V-Ar/72 e**)eE 

Take the transition 

(u,v) -» (u'.v' + V'). 

Then the equations (|5| are rewritten as (drop the primes) 

dip 



dt 

On 
dt 



[iiAip + Xiip + Xiip* - 21/44/ — -u - 27; 



72 



^ l^l 2 - 2v A uip - 7 2 |V ; | 2 V', 



72 



-/i 2 A^u + A 2 Au + 2i/ 4 4 / + A[(-i/ 2 + 3i/ 3 JST)ir + i/&u d + v 4 \ 

V 72 



udx = 0, 



where -0 = -01 + ^2 



We consider the transition of (22 1 beyond E. The linear operator of (22 1 is given by 

/ 



2Ai -2 Vi « f=f± \ 
V 72 



B = Ai 

y2i/ 4v /^A A(- M2 A + A 2 )y 
Restricted to its first eigenspace 

El = [(V'l, "02, U) = (j/1,2/2,2/3) COS ^ I (2/1,2/2,2/3) G R 3 } , 



(21) 



(22) 



(23) 
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the linear B is given by 



D\ 



2Ai 




-2v A 
Ai 



72 



72 L 2 \ L 2 



The eigenvalues f3\, (3 2 and {3$ of B\e 1 are given by 



ft 



Ai 



2L 2 



I/ 2 



+ A : 



± 



Ai 



7T Z / /X 2 7T^ 



2L 2 V L 2 



+ A 2 



2Al7T 2 / /i 2 7T 
"L 2 " IT 2 " 



+ A 2 



ft = Xi < 0. 
We know that 

By d9l, we have 



ft < 0, ft < ft VT < T c: 



A 2 = *l - 2^ 2 X + 3/y 3 ^ 2 = 0i (T + <r ) - 2v 2 (l - X)X 
2v 2 



= 9 1 [T + a a - -^X(l-X) 



h (T - T c2 ) - 



L 2 



Hence the transition curve f3\ = is given by 



which is equivalent to 



L 2 



+ A 2 



?i (r-T c2 ) 



2^| 

72 



2^ 

72 



Hence the transition curve /?i = is given by 

ft = T = T c2 



0i 7 2 i? 6>iiv 6*172 



2^f 

72 



(24) 



(25) 



Then L\ defined by (16) is the critical length scale to make the transition curve fl\ = achieving its maximum at 
(T,X) = (0,1/2). 

By definition, it is then easy to see that if L < L\, no phase separation occurs at any temperature although phase 
transition for He does occur as the temperature decreases below certain critical temperature. 

Step 3. Now we calculate the length scale L 2 where the critical curves ft = and ft = are tangent to each 
other, i.e., they interact at exactly one point. By definition, we need to solve L such that the equation 

_ (1 _ X) = _ X(1 _ X) _ C70 ___ 



has exactly one solution. Hence, it is easy to derive that the formula L 2 is as given by (17 1. 
By comparing L\ and L 2 , we have 

L 2 > L\ if and only if < —. 

62 R 



(26) 
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Now we the length scale L3 where the critical curves Pi = and fii = are tangent to each other, i.e., we need to 
solve L such that the equation 



ft 1 -* 



2 a 



X(l-X)-* 



has exactly one solution. The formula L3 is as given by (19). 

Step 4. Now we return to prove Assertions (l)-(4). First, we consider the case where ^ < j[, i- e - 5 L 2 > L\. In 
this case, for L\ < L < L2, second transition curve /3i = is as shown in Figure [3) Namely, as one decreases the 
temperature T entering from region D2 into -D3 through the curve (3\ — 0, phase separation occurs. 

By (251, solving 



-X{l-X)-o,-— 2 



2vl 
6*172 



= 



gives the least and biggest mol fractions X\ < X2 defined by (19 1, and Assertion (2) follows. Other assertions can be 
proved in the same fashion. 

The proof of the theorem is complete. ■ 



IV. PHYSICAL CONCLUSIONS 



We have introduced a dynamical Ginzburg-Landau phase transition/separation model for the mixture of liquid 
helium-3 and helium-4. In this model, we use an order parameter i/j for the phase transition of liquid 4 He between the 
normal and superfluid states, and the mol fraction u for liquid 3 He. As u is a conserved quantity, a Cahn-Hilliard type 
equation is needed for u, and a Ginzburg-Landau type equation is needed for the order parameter. The interactions of 
this quantities arc built into the system naturally by using a unified dynamical Ginzburg-Landau model for equilibrium 
phase transitions, where the dynamic model is derived as a a gradient-type flow as outlined in the appendix. 

This analysis of the model established enables us to give a detailed study on the A-phase transition and the phase 
separation between liquid 3 He and 4 He. In particular, we derived three critical length scales L\ < L2 < L3 with the 
following conclusions: 

1) For L < Li, there is only A-phase transition for 4 He and no phase separation between 3 He and 4 He, as shown 
in Figure [2j a). 

2) For L\ < L < L2, there is no triple points, and phase separation occurs as a second phase transition after the 
A-transition when the mol fraction is between two critical values, as shown in Figure [3] 

3) For L2 < L < L3, the A-transition is always the first transition, as shown in Figure |4j 

4) For L3 < L, both the A-transition and the phase separation can appear as either the first transition or the 
second transition depending on the mol fractions, as shown in Figure [5] In this case, when the phase separation 
is the first transition, the separation mechanism is the same as a typical binary system as described in great 
detail by the authors in pQ. 

Also, the A-transition is always second-order. The derive theoretical phase diagram based on our analysis agrees 
with classical phase diagram, and it is hoped that the study here will lead to a better understanding of mature of 
superfluids. Finally, we remark that the order of second transition is mathematically more challenging, and will be 
reported elsewhere. 



APPENDIX A: DYNAMIC GINZBURG-LANDAU MODELS FOR EQUILIBRIUM PHASE 

TRANSITIONS 

In this section, we recall a unified time-dependent Ginzburg-Landau theory for modeling equilibrium phase transi- 
tions in statistical physics; see also [5J[B]. 

Consider a thermal system with a control parameter A. By the mathematical characterization of gradient systems 
and the le Chatelier principle, for a system with thermodynamic potential H(u,\), the governing equations are 
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essentially determined by the functional Ti(u, A). When the order parameters (iti, • • • , u m ) are nonconserved variables, 
i.e., the integers 

Ui(x,t)dx = a,i(t) ^ constant. 



then the time-dependent equations are given by 



duj 
dt 

for 1 < i < m, where fli > and $i satisfy 



-A— W(<u,A) + $;(u,V«,A), 

OUi 



f 6 

/ V H(u,\)dx = 0. 
Jn om 



(Al) 



(A2) 



The condition ( A2 ) is required by the Le Chatelier principle. In the concrete problem, the terms <&i can be determined 
by physical laws and (A2|. We remark here that following the le Chatelier principle, one should have an inequality 
constraint. However physical systems often obey most simplified rules, as many existing models for specific problems 
are consistent with the equality constraint here. This remark applies to the constraint (A8l below as well. 

When the order parameters are the number density and the system has no material exchange with the external, 
then Uj (1 < j < m) are conserved, i.e., 



/ Uj (x, t) dx = constant. 
Jn 



This conservation law requires a continuity equation 

duj_ = 
dt 

where Jj(u, A) is the flux of component Uj, satisfying 



-V- Jj(u,X), 



where fii is the chemical potential of component it;, 



SUj 



H(u,X) — 4>j(u, Vu, A) , 



(A3) 



(A4) 



(A5) 



(A6) 



and <j)j(u,\) is a function depending on the other components u$ (i =/= j). Thus, from (A4l-(A6l we obtain the 
dynamical equations as follows 



dt 



-7i(u, A) — 4>j (u, Vu, A) 



for 1 < j < m, where (3j > are constants, and 4>j satisfy 

/ n £^r^A )ffa = o. 



(A7) 



(A8) 



When m = 1, i.e., the system is a binary system, consisting of two components A and B, then the term <f>j = 0. 
The above model covers the classical Cahn-Hilliard model. It is worth mentioning that for multi-component systems, 
these <fij play an important rule in deriving good time-dependent models. 

If the order parameters (iti, • • • , rife) are coupled to the conserved variables (ufe+i, • • • ,u m ), then the dynamical 
equations are 



duj 
dt 

dt 



-f3 i —H(u,\) + <P l ( u yu,\), 
oui 



(3 3 A 



5u 



Ti.(u, A) — 4>j(u, Vu, A) 



(A9) 
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for 1 < i < k and k + 1 < j < m. Here <I>i and <j>j satisfy (A2| and (A8l, respectively. 

The model ( A9 1 we derive here gives a general form of the governing equations to thermodynamic phase transitions, 
and will play crucial role in studying the dynamics of equilibrium phase transitions in statistical physics. 
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